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Boundary conditions changing operators have played an important role in conformal 
field theory. Here, we study their equivalent in the case where a mass scale is introduced, 
in an integrable way, either in the bulk or at the boundary. More precisely, we propose 
an axiomatic approach to determine the general scalar products . . . , 9 m \\9[, . . . , 0' n ) a 
between asymptotic states in the Hilbert spaces with a and b boundary conditions re- 
spectively, and compute these scalar products explicitely in the case of the Ising and 
sinh-Gordon models with a mass and a boundary interaction. These quantities can be 
used to study statistical systems with inhomogeneous boundary conditions, and, more in- 
terestingly maybe, dynamical problems in quantum impurity problems. As an example, 
we obtain a series of new exact results for the transition probability in the double well 
problem of dissipative quantum mechanics. 



01/98 



Packard Fellow 



1. Introduction 

The concept of boundary conditions changing operators has been of crucial importance 
in the analysis of boundary conformal field theories, as well as in their applications to 
quantum impurity problems. In the last few years, important progress has been made in 
extending the solution of conformal field theories to theories perturbed either by a bulk or 
by a boundary operator. The latter have extremely interesting applications to the study 
of flows in quantum impurity problems. 

As we will show here, it is possible to introduce boundary conditions changing opera- 
tors even in the case where there is a mass scale, either in the bulk or at the boundary. It 
is not completely clear what the formal use of these objects (eg as illustrated in 0] for the 
conformal case) might be, but they certainly do have applications. From a 2D, statistical 
mechanics, point of view, one might wonder what is the effect of having different parts of 
the boundary with different boundary conditions 0, describing, for instance a situation 
with a (classical) "impurity" on the boundary. From the 1+1 point of view of quantum 
impurity problems, one might wish to describe situations where the coupling to the im- 
purity is changed at some particular time, and one is interested in the subsequent time 
evolution of the degrees of freedom. In fact, one of the key observables in the two state 
problem of dissipative quantum mechanics, the well known quantity "P(i)" (see below), is 
essentially defined in that fashion [[|. 

In the conformal situation, only a discrete set of conformal boundary conditions are 
available, and the effect of switching from one to the other is described by the insertion 
on the boundary of the appropriate operator. For instance, in the Ising model, one goes 
from free to fixed spins by inserting the conformal operator $12 of weight h = and 
from fixed up to fixed down by inserting $13 of weight h = \ |IJ (where the $ rs are 
the usual degenerate conformal fields ||). We are interested in the more general case 
where conformal invariance is broken at the boundary, and also maybe in the bulk. As 
an example, we can consider a situation where, in the Ising model, we apply a boundary 
magnetic field h a for y < 0, hb for y > (we use coordinates x and y to describe the two 
dimensional space and the boundary sits at x = in this paper), and in addition, T ^ T c 
in the bulk. In general, one does not expect this situation to be described by the insertion 
of a simple operator at y = anymore. Nevertheless, the change of boundary conditions 
can be fully characterized in the factorized scattering description (the Ising model with 
a boundary field being integrable) by scalar products of asymptotic states in the Hilbert 
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space with h a and h\> (see below) , and the "operator" inserted at y = can thus be written, 
in principle, in terms of the Faddeev Zamolodchikov algebra. The problem then reduces to 
the determination of scalar products (we also call them transition factors), and is similar 
in nature to the problem of determining form- factors in integrable theories f|. 

In this paper, we mostly restrict to the Ising case, where the computations are already 
moderately complicated, and obtain the complete solution to the problem with changing 
boundary magnetic fields in the massive theory. In section 2, we discuss the general 
problem of form-factors in the crossed channel. We determine the scalar products of 
interest for h a hb > in section 3. In section 4, we discuss the limit where the bulk is 
massless, including the conformal invariant case. In section 5, we discuss the case h a hb < 
which require a different treatment. Section 6 is devoted to extending some of these results 
to the sinh-Gordon model, where the non trivial bulk scattering matrix introduces further 
complications. In section 7, based on some reasonable conjectures, we apply and generalize 
our results to the determination of P(t) in dissipative quantum mechanics, for arbitrary 
value of the dissipation parameter g. 

Some of the results presented in this paper have appeared in a shorter version ||. 



2. Form factors in the crossed channel 

Consider the massive Ising model defined in the half plane x G (— oo, 0], y G (—00,00), 
in the presence of a boundary magnetic field. The action reads 

/0 pod 2 r°° /"oo 

dx / dy apF(x,y) + — / dy [(V'V') {x = 0) + ad] + h / dyas{y)- (2.1) 
-00 J —00 * J — 00 J— 00 

Here apF is the usual massive free Majorana fermion action, a is a boundary fermion 
satisfying a 2 = 1, is the boundary spin operator, which coincides with | (ip + -0) a. 

As discussed in the problem can be studied from the point of view of the direct 
channel with x taken as time and y as space. In that case the Hilbert space is the usual 
one, and the boundary represented by a boundary state. On the other hand, in the crossed 
channel, one has a new Hilbert space for the theory on the half line. Setting z = x + iy, 
the one point function of the energy is easy to compute in the direct channel. Using the 
boundary state formula 



\B) = exp 



00 jg 

—K(9)Z*(-9)Z*(9) 
-00 47r 



|0>, (2.2) 
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where K is related to the reflection matrix through K{9) = R ( ^ — 9), and the two particle 
(bulk) form factorB 

(0\e\6i, 9 2 ) = im sinh {^-^j , (2.3) 

one obtainsl 

(e(z,z)) = -i— [ —K(9)smh9e 2mxcoshd . (2.4) 
2 J 2ir 

An expression in the crossed channel follows by moving the contour of integration in the 
variable 9 (w = iz) 

(e(w,w)} = - [ —R(9)cosh9e- 2imxsinhe . (2.5) 
2 J 2% 

There is no simple way to compute this expression directly in the crossed channel, except 
by solving the problem explicitely, writing mode expansions for the fermions (see appendix 
A). In particular, ( |2.5| ), which can be considered as a zero particle form-factor, does not, 
as far as we know, follow from any known nice form of the ground state of the theory with 
boundary. 

Using the mode decomposition of fermions, it is possible to construct form factors 
with more particles. Because of the boundary, the states at 9 and —9 are related: one has 
Z*(9)\0) a = R(9)Z*(—9)\0) a (here |0) a denotes the ground state with boundary condition 
of type a at x = 0, normalized such that a (0|0) a = 1. We drop the label a when it is not 
necessary). We denote by ||0) o the one particle asymptotic state on the half line Z*(9)\0) a 
i. Then one has 



\e(w,w)\\9i92) = im e~ 



iy (cosh Q\ +cosh 02 ) 



Yl [1 + R(9 k )t k ] Sinh ( — J ^(smh^+sinh^)* 

fc=l,2 V 2 / 



(2.6) 



where t% is an operator acting on functions of many variables as follows 

Ug{---9i---)=g{ 0<-..), t 2 = l, t i t j =t j t i . (2.7) 

This expression, call it J- a a(0i&2), obeys the relations 

Faa{9l, 9 2 ) — R{9\)J : aa {—9\, 9 2 ) = R(92)J I 'aa(9l, — 6> 2 ) = —J r aa{9 2 , #1, 



1 Here, our normalization for bulk asymptotic states is is (#i|#2) = 2tt8{0\ — 62). 

2 Rapidity integrals, unless specified, run from —00 to 00. 

3 The normalization is the same as in the bulk: (9i\\02) = 2tx5{6\ — 62). 
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This will allow the use of integrals running from — oo to oo in the computation of correlators. 

Observe that ( |2.6| ) does not have a pole at 0% = 0\ + in, as is only natural for the 
energy operator. Therefore, in the crossed channel, there is no relation between the two 
particle and the zero particle form-factor of the energy via kinematic poles. 

In order to better understand and generalize these form factors, let us follow an 
axiomatic approach to try to compute them. Recently, axioms for form factors in the 
cross-channel in the presence of a boundary have been written for the XXZ spin chain || . 
For notations sake, we define the form factor in the cross-channel with boundary conditions 
of type a to be 

F aa {9i, -, On) = a (0\OZ*(9 1 ), Z*(0 n )\0) a . (2.9) 

Following the authors of || and taking the continuum limit of their axioms we find, in our 
specific case that 

Faa{9li ■ ■ ■ , Oi, ' ' ' n ) = S{9i — i + \) T aa {0 1 , • ■ • , Oi, ■ ■ ■ n ) (2-10) 

Faa{0l, ■ ■ ■ , Oi + 2ni, ■ • • n ) = T aa \0\i • • • , Oi, - • • O n )S(0i-l — Oi) ■ ■ ■ S[0l — Oi) 

x Rl{0i + i^Si-Oi -Oi)--- Si-di-t - 0i)S(-0 i+1 -Oi)--- S(-6 n - Oi) (2.11) 

xR* a (0 i )S(0 n -0i)---S(0i +1 -0i) 

The kinematic pole equation gives the residue of the form-factor as 0j = 0± + in, and it is 
given by 

ReS.F aa (6>i, ■ • - , Oi, • ■ ■ , n ) = lT a a{0\ ■ ■ - Oj ■ ■ ■ n ) 

[S(#! - Oj) ■ ■ ■ Sip,-! - Oj) - S(0j - J+1 ) ■ ■ ■ S(0j - n )R a (0j) (2.12) 
xS(0 n + 0j)...S(0 2 + 0j)R* a (0j)}. 

Other residues are also generated from the bound states of the S and R matrices, if any. 

We have written the axioms without internal indices, which are implicit in the previous 
equations. In the case of diagonal scattering, the R matrix actually drops out from these 
equations due to RR* = 1. The only dependence appears in the last, reflection equation 

^aa(0l, ■ ■ ■ , Oi, - ■ ■ , n ) =S(0i — ■ ■ ■ S[0i — n )R a (0i) S (—0 n — Oi) 

(2.13) 

X • • • S{— Qi+l — O^TaaiOl, • • • , —Oi, ■ ■ ■ , n )- 

Up until now, solutions to these axioms have not been written and this is the task we plan 
to undertake here. In fact, we will first generalise these form factors slightly, and this will 
lead to a more general set of axioms. 
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3. The case of mixed boundary conditions 

3.1. The two particle case 

We now consider the case of inhomogeneous boundary: we suppose that the boundary 
conditions for y G (—00, 0] are described by a boundary magnetic field h a , and for y e [0, 00) 
by a field h b . Suppose we wish to compute again (e(w,w)) in the cross-channel (in the 
direct channel, no expression for the boundary state \B) is available). To do so, we have 
to be careful that the asymptotic states will depend on the boundary conditions in a more 
intricate way than through the simple relation \\6) a = R a (6)\ \ — 9) a . Indeed, the scalar 
products 

o\ v nj ■ ■ ■ i Ol\\0n+1 ■ ■ ■ 9 n -\-m/a 

have to be non zero in general, even for disjoint sets of rapidities (m even), since the 
asymptotic states provide a complete set of states for any given boundary condition. As a 
result, the one-point function of the energy with mixed boundary conditions will read (for 
y>0) 

/•oo in jq 

(e(w,w)) ba = / ^ T 1 {6(O||e(«;,«;)||^i02) b 6(Mi||O) a 

Jo 8tt- (3.1) 

+ b (0\\e(w,w)\\0) b 6 <0||0> a }/ 6 (0||0> a . 
Let us now introduce the quantities 

b(O20l\\O}a 

6<0||0) a 

and 



G(9 1 ,e 2 ), (3.2) 



b (0\\0) a ~ ( ' ^' 

To determine these scalar products we sometimes call transition factors, we use, as in the 
well known case of operators in the bulk, an axiomatic approach || . It is then sraightfor- 
ward to write the first set of axioms 

G(d 1 ,e 2 ) = -G(e 2l 6 1 ) = R* b (e 1 )G(-e ll e 2 ) = R* b (e 2 )G(e ll ~e 2 ). (3.4) 

To complement these, we need a residue condition. It is easily obtained as follows. Consider 
the one point function of the energy. It does have different expressions for y > 0, where 

(e) ha = ^j ^R b {0)^ee-^^-im J ^^(^)sinh^ 
G(9i, 9 2 ) exp[zm(— sinh6*i + sinh^)^] exp[— m(cosh6'i + cosh 9 2 )y] 
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and for y < 0, where 

2, J 111 J 07T Z Z 

F(d±,d2) exp[*m(sinh0i — sinh 6*2)0;] exp[m(cosh6>i + cosh 6*2)2/], 



where we used the symmetry relations ( |3.4j) to write all integrals on the interval (—00, 00 
By requiring that the second expression is the analytic continuation of the first, we find, 
moving the contours of integration, that G must have a simple pole for 62 = 0\ — in, with 
residue as 6*2 — > #i — ztt, 

Res G(9, 6-m) = -i (l - tjM} . (3.5) 

This is equivalent (see below) to having a simple pole as 6*2 — > 6>i + Z7r 

Res G(0, + ztt) = -i (l - |^|^) . (3.6) 

Similarly, F has poles for the same values 6*2 = 0\ ± z-zr, with residues switched. 
Recall that the -R matrix reads 

/ . 7T 6\ n a — i sinh 6* 

it a = 2 tanh — -, (3.7) 

V 4 2y K a + ^sinh6 !, v ; 

with n a = 1 — /i^/(47r). It obeys the relations 

R(e + i7T) = - 1 ^-r = -R*(e). (3.8) 

In the following, we parametrize K a = — cosh^a (and similarly for Kb)- It is useful to 
observe that 

— — = 1 tanh 1 — % tanh 1— , (o.9) 

Ka + zsinhfl V 2 4/ V 2 4/' v; 

Let us now try to determine the functions F, G. To do so, consider the equation 

*(6)*(8 + in) = 1 _ e (3.10) 

cosh - if) cosh (z-f* + ^) 



4 When doing so, one has to use the normalization (#1 1 1 #2) = 2n[S(9i —62)+ R*{6i)8{6\ + 62)]- 
The presence of two delta terms is compensated by factors | one has to introduce to get the in- 
tegrations running over the whole real axis. 
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Its simplest solution is 



n 



r(t + n-<^) 



cos h(^-f) n y o r(! + --^) 



r(f + n + <^)r 



n 



(3.11) 



r(! + - + ^)r(f + n-^)r(| + n + ^) 

or, using an integral representation, 



cosh - f ) 



exp 



oo 7. j Q-Qo- t i b t 

dte 1 £ — e* 2^ 1 
.oo t 4 cosh I sinh I 



(3.12) 



The function $ has no poles in the physical strip Im 9 £ [0, 7r], and a single pole 6> = 9 b — 1 f 
in Im 6> £ [— 7r, 0]. It satisfies the following identities 



**(0|06A) = *(-0|-06,-0a) 



together with 



$(0 + i7r|0 6 ,0 a ) = 



1 



tanh(^ -if) 

$(0|0 O ,0&)- 



(3.13) 



(3.14) 



itanh(^^ - if) 

The asymptotics of $ can be worked out easily from the integral representation. If 9 — 
9 a ,b —> oo, one has $ m 2cue~ 6 '/ 2 e^ a+6 ' b ^ 4 (a; = e l7F//4 ), while in the opposite case 9 — 9 a ,b —> 
-oo, one has $ « 2cu~V /2 e~ (ea+eb)/4 . 
We now introduce the quantity 



f(9\e b , 9 a ) = y/-i{K a -K b )${9\9 b , 9 a )$(9\ - 9 b , -6 a ). 



(3.15) 



It is important to stress that, despite the notation, this function is independently an even 
function of 9 a and 9 b - By convention, we always chose these variables to have positive real 
part in what follows. The function / obeys 



f(9)f(9 + i7T)=- 



(«a - «&) 



icosh - if) cosh - if) cosh + if) cosh + f ) 
sinh# I ' 



together with 



f{e + m\e b ,6 a ) = f*{6\6 a ,0 b ), 



(3.16) 
(3.17) 
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and 



f(e + 2in\e b ,e a )=R a R* b (9)f(e\9 b ,9 a ), 
n b — i sinh 9 



f(0\9 b ,9 a ) 



-J(-9\9 b ,9 a ) 



(3.18) 



(3.19) 



K b + i sinh 9 ' 

Based on these identities, we can now write the minimal solution for the form factor 
axioms fltP] ) and (|3T6|), 

n b + i sinh 9i 1 



G(9 1 ,9 2 ) =- l[ f(9i) sinhfli 



i=l,2 



K b -isinh9i C osh(%+if) 



x tanh ( — !— : — -| tanh ' 



(3.20) 



2 J \ 2 

From the expression of G together with ( |3.14j ), one easily checks that 

G(9 1 + 2m, 9 2 \9 b , 9 a ) = -R a (9 1 )R* b (9 1 )G(9 1 , 9 2 \9 b , 9 a ), 



(3.21) 



so the residue axiom ( |3.6| ) is also satisfied. The other form-factor F can also be easily 
obtained. For this, observe that 

F{9i,9 2 \9 a ,9 b ) = [G(9i,9 2 \9 b ,9 a )]* 

Using (J3.8[),(|3TTBD, the crossing identity follows: 



that is, explicitely 



F(9i, 9 2 ) — —G{6\ — m, 9 2 — in), 



b(0\e 1 ,9 2 ) a = b (9 1 -m,9 2 -in\0) a . 



(3.22) 



It follows for instance that 

F(9 1 + 2in, 9 2 \9 b , 9 a ) = -R* a {9 l )R b {9 1 )G{9 l , 9 2 \9 b , 9 a ). 

Finally, from ( |3.22| ), we deduce the result 

b{9 2 \9\)c 



,<0|0) a 



H(9 1 ,9 2 )=G(9 1 -in,9 2 ) 



(3.23) 



(3.24) 



It follows that H has a pole as 9 2 — ► 9\, with residue — i [1 — R a Rl(9\)]. It also has a pole 
as 9 2 — > —9\ with residue i [R b (9i) — R a (0i)]. As in the bulk case, these poles indicate 
that ( p.24| ) holds only for "distinct" rapidities ||; the term 



<S(0i - 9 2 ) [1 + R a Rt{9 x )] + 5(9, + 9 2 ) [RM + R b {9{)] , 



(3.25) 
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has to be added at coincident rapidities. This quantity, H, is probably the easiest to 
understand intuitively. For instance, compare our results with elementary computations 
for bare particles. With boundary condition a at x = (and a trivial boundary condition 
at x = —L), a state \9 a ) is associated a one particle wave function 

ip a (x, 6) = -= ( e imsinhex + R a (6)e- irns ' inhdx ) . 
v2 

The quantization condition is e 2imsinheL R a (9) = 1, and the normalization 

r° i 

/ \iP\ 2 dx = L — TmRJd) ^L,L» 1. 

J_ L msmhO 

One then finds the scalar product 

'Rt(e 2 )R a (e 1 )-i R* b (e 2 ) - Ratfi 



2m 



sinh 9 2 — sinh 9\ sinh 9 2 + sinh 9 



where we used the quantization conditions. This expression has the same poles, with the 
same residues, as in ( |3.24| ). The complication for the physical theory is of course due to 
the fact that the ground state shifts with different boundary conditions, hence "dressing" 
the previous expression. 

Inspection shows that G has a simple pole at 9i — ^ and at 0j = ±6>b + (the pole 
of $ at 9 = 9b — is cancelled by the Kb — i sinh 9 term in G) . Similarly, F has a simple 
pole at 9i = — ^ and at $i = ±9 a — -y. The poles of G are thus the same as the poles of 
Rl, the poles of F the same as those of R a . 

Notice that the residue condition does not prevent the quantities F, G from having an 
added part without kinematic pole, with an expression similar to the ones written so far, 
but with the tanh of the sum and difference of rapidities replaced by the sinh function. 
However, dimensional analysis, together with the requirement of having a well defined 
massless limit and comparison with perturbation theory, exclude such a term. 

Notice also that, like the R matrices, the scalar products depend only on squares of 
magnetic fields. It is easy to check that this is true to the first orders in perturbation 
theory, see the appendix. It is important to stress that our results apply in fact to the 
case of fields of the same sign; when h a hb < 0, slightly different formulas have to be used, 
which we discuss below. 
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We can finally compare some of these results with perturbation theory. Consider the 
situation where the two boundary fields are very small; at leading order, the R matrices and 
the functions $ are evaluated for vanishing fields, ie a = 9 b = itt. One has in particular 

1 



$((9|z7r, in) 



cosh(f +tf)" 

Putting these values back in the formula, we find, after some algebra 

G(6>i, 62) oc tanh — - tanh 1 — - cosh ( + i — 

1^2 tt\ sinh 6*i sinh 6*2 
x cosh — + 1 — 



2 ~4J cosh 2 6 X cosh 2 6 2 ' 

On the other hand, the scalar product can be computed from straightforward perturbation 
theory. The fermion propagators are obtained by solving the equations of motion, using 
formula (4.3) of 0, which provide directly integral representations. After some lengthy 
algebra, one finds a result in agreement with the previous formula for G. 

3.2. The general case 

Following the algebraic approach of [@ we see that we could as well put different 
boundary conditions in the left vacuum to that of the right vacuum in the form factors 

FbaWi, ■■■,e n )= b (O|0Z*(0i) • ■■Z*(d n )\0) a (3.26) 

Then the residue axiom for the two particle form factor becomes 

Res^ ba (d 1 ,e 2 )\0 2=ei+in = i[i-R a (e 1 )Rue 1 )] b (o\) a (3.27) 

and can be generalised to the n particles form factor in the obvious fashion. Here it is 
important to emphasize that the .F's are not form factors in a theory with inhomogeneous 
boundary conditions, but rather "transition form factors" between two theories having 
different boundary conditions (and thus different Hilbert spaces) . These axioms are general 
for operators which are local with respect to the Fadeev-Zamolodchikov operators. In 
particular, if we choose the case O = 1, ie the identity, then the previous axioms will yield 
the scalar products (or transition factors) between asymptotic states in Hilbert spaces with 
different boundary conditions. The other axioms have to be modified in a similar way; for 
instance, one has 

T ha {Qx + 2m, e 2 ) = s(e 2 - e 1 )R* a (e 1 )S(-e 1 - e 2 )R* b (9 l + m)T ba {e t , e 2 ). (3.28) 
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The modification of the exchange type reiations is trivial. 

In the special case of the Ising model where S = — 1, we get the simpler equations (for 
a local operator) 



T ha {- • • 0< + 2wr • • •) = -R*(ei)R b (0i)Fba(- ■ ■ &i ■ ■ ■) 
Rese^+inFbaiOt, ■ ■ ■ , d t , ■ ■ ■ 9 n ) = i(-iy +1 [1 - R a (9 1 )R* b {9 1 )] 

X J~ba\6l-> •••%)•••) On) 
Fba{- ■ ■ Oi- ■ •) = R a {9i)Fba{- ■ ■ — 9{ ■ ■ ■). 



(3.29) 



Given the function / constructed before, it is now easy to write the generalisation of 

f>(02n, • • -0l|O) a T-r . , .Kb + ISwhdi 1 

— — =c n \\Bmhdif{O i ) —r-TT, r JW- , ^ 

6(0|0) a ^ k & - zsmh^ cosh + l -f) 

9 - — 9 ■ 9 + 9 ■ (3.30) 

x I [ tanh — tanh — - 

11 2 2 

with c n a proper normalisation determined by the residue condition, c n = (f)™. The 
various other components are obtained by crossing, in particular 

&(O|0i,---02n}o / a„ tt ■ , a . ,K b -isinh9i 1 

— — = -1 c n \\s\nh9 l f{9 i -m) — -jj. ^ 

&(0|0) a - L .- L «t + zsinh0i cosh - ^) 

Q. — Q . 9+9- (3.31) 

x TT tanh — tanh — - 

11 2 2 

These formulas also agree with the first non trivial order in perturbation theory. 

3.3. Compatibility of form-factors and transition factors 

It must be clear that although asymptotic states with different boundary conditions 
have some non trivial scalar products, the form-factors of physical operators have the same 
general form with any boundary conditions, ie all they see from the asymptotic states is 
their property \\9 >= R(8)\ \ — 9 >. This is made possible by the axioms satisfied by the 
transition factors, in particular the pole condition. Let us illustrate this in the case of the 
fermion operator. We have 

a(0W\9)a = oo (e e ' 2 + R a {9)e~ e ' 2 ) . (3.32) 

On the other hand, we can compute this form-factor by inserting a complete set of states 
with b boundary conditions on the left and on the right. Using the behaviour of transition 
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factors under reflections, this allows us to reexpress (|3.32j ) as a sum of two types of terms. 
The first type is 



UJ 



+ 



^V^ 2 [ a <o|ou^> a 



d9 2 d9 3 
(4tt) 2 2! 



(3.33) 



and the second type is 



UJ 



+ 



jgj e fli/2 

2n 

d9 2 d9 3 d9 4 



a (0\9 2 9 s ) bb (9 3 9 2 9 1 \9) a + 



df) 

^ a (o\6ie 2 ) bb (d 2 \d) a 

An 



a (0\6 1 6 2 9 3 9 4 ) bb (6 4 9 3 6 2 \6) a + 



(3.34) 



(4tt) 3 3! 

In all these formulas integrals are regulated by taking principal parts. Let us now consider 
the second sum ( |3.34| ). The 9\ integral has no singularity on the real axis; we move the 
contour to Im (#1) = 7r; by doing so we do not encounter any singularity since the transition 
factors with 9\ on the right have poles only in the lower half plane. By using crossing and 
the completude relation with b boundary condition, this gives rise to 10 (e e / 2 + R a (9)e~ 6 / 2 ) , 
up to two corrections. First, the correct crossing formula is a (0|6>i + in, 9 2 ) b = a (9i\9 2 )f c , 
where on the right we mean the disconnected transition factor, without the added delta 
function part at coincident rapidities. On the other hand, the completude relation works 
with a (9i\9 2 ) b : this means that we have to subtract a first correction, which reads using 
( |3.25| ) , together with a bit of algebra, like (|3.33| ) but with the integrand for the 6>i variable 
multiplied by 1 + \^{9i) + §^ Second, after crossing, we get a 9\ integral that 

runs over the whole real axis, avoiding the poles at 9 2 and —9 2 by going under them. 
The completude relation on the other hand requires the principal part integrals, so we 
have to subtract the second correction, given by the pole contributions. This gives rise 
again to a formula like ( |3.33| ) but with the integrand for the 9\ variable multiplied by 



1 _ l3a(Q 

1 2 R h \ u 



\-r-{9i)i where we used the residue conditions for H. Adding up the two 



corrections, we see that the second type of terms (|3.34| ) is nothing but (|3.32|) minus the 
first type of terms (|3.33|) . Hence, adding up ( p.33|) and ( |3.34j ) gives (|3.32|) , and we recover 
our form factor. 



4. Massless limit and boundary conditions changing operators 

We now consider the case where the bulk is in fact massless. This is conveniently 



adressed within the framework of massless scattering , |Kj , which is formally obtained 
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by letting the physical mass m — > 0, together with the rapidities 9 — > ±00. We set 
m = % e _e> (/jl an arbitrary scale, taken equal to unity in what follows), 9 = ±0 ± f3, 
to parametrize energy and momentum e = ±p = //e^ 3 . The ± sign corresponds to R (L) 
movers. 

We first consider time propagation in the x-direction (direct channel). Then, L and 
R movers are defined by the conventions, eg for the form-factors of the fermions 

(0| Vi?(x, y) \(3) R = ue^ 2 exp [e?(x + iy)] 

(4.1) 

(0\Mx, v)\P)l = we^ 2 exp [ e P(x - iy)} . 

Consider then the one point function of the energy operator, e = itpL^R- Using the 
expression for the boundary state 

i s > = E^ / U^ K ^-^)zm)zm)\o), (4.2) 

n=0 ' ^ i 

one finds (here, z = x + iy) 

(e(z,z)) = (0\e(z,z)\B) = -i J ^-K((3 - (3 B )e^e 2e ^. (4.3) 

In the case of fixed boundary conditions K = i, and one finds (e) = — 4^:. For free 
boundary conditions, K = —i, and one finds the opposite result. 



In contrast with the massive case, where the expression (O) could not be obtained 
directly in the cross-channel, such a computation is now possible. This is because we can 
represent asymptotic states as superpositions of left and right moving parts, the simplest 
ones beingS 

\\(3)=\(3) R + R(f3)\(3) L 

WW =\Pl,fc)RR + R(p2)\PMRL + R(Pl)\PMLR (4.4) 
+ R(Pl)R{02)\0102)LL, 

and use the LR factorized form of conformal operators. Here, the meaning of L and R is 
encoded in the new dependence of form-factors 

(0| \ip R (x, y) \\P)r = ue^ 2 exp [e^(-y + ix)] 

(4.5) 

<0| \^ L (x, V) \\P)l = ^' 2 exp [ e P(-y - ix)] . 



These states are normalized such that, for instance, (/3i||/32) = 2ir5(/3i — P2) 
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Using the general crossing formula written earlier, the one point function of the energy 
reads then (w = iz) 



(e(to,tD)> = -i / ^(OlI^^H^^HVflHIlO) = / - p B )ePe-»** . (4.6) 



A rotation of the countour shows that the expressions (|4.3| ) and ( |4.6|) are identical, with 
the identification i*T(/3) = -R(^f- — ft), as in the massive case. The massless reflection 
matrices actually follow from the massive ones by setting 9 b ~ + they are of the 
form = i tanh (S^k - f ) . 

We will define massless form-factors as form-factors on the previous asymptotic states. 
They follow from the massive case by taking the infinite rapidity limit 9 = ±0 ± f3 for 
each rapidity. For instance 



(0| \e{w, w) | |/3i/3 2 ) = ie pl/2 e p2/2 {R((3 2 ) exp[e pl {-y + ix) + e p2 (-y - 



IX) 



(4.7) 

R(/3i) exp[e pl (-y - ix) + e p2 (-y + ix)] } . 

The scalar products between asymptotic states with different boundary conditions follow 
by the same limiting procedure. Alternatively, one could obtain them by writing a similar 
set of massless axioms. As an example, in the massless case, the expressions for the one 
point function of the energy are 



G([3 u f3 2 )e pl/2 e p2/2 exp [e pl (-y - ix) + e p2 (-y + i 



(4.8) 



for y > 0,and 



(e(x, y )) ba = J^-R a ((3)e p e-^ +i J ^KVh 
G(f3 u p 2 )e pl / 2 e p2 / 2 exp [e pl (y + ix) + e p2 (y - i 



(4.9) 



for y < (where we used the antisymmetry of F, G) leading to the same condition as 
with now the massless reflection matrices and in the variable (3. Using that f(0\9 b , 9 a 



2e~ 2 ^Wsinh pb pa $((3\/3 b , (3 a ) gives the result 



G(/3 1 /3 2 ) = -^sinh^ y ^ 



n m\hM^w 



tanh- — ^ (4.10) 
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where we have set T a (&) = e^ a{h) . 

We now consider the case where the mixed boundary conditions are conformal invari- 
ant: we chose free boundary conditions (F) for y < and fixed boundary conditions (+) 
for y > 0. This is described in the previous formalism by taking the limits /3 a — > — oo and 
[3b — > oo. From previous formulas, we find 

+ (0\0)f =' taDh — 2~ • ( 4n ) 
The one point function of the energy reads then 

(e(w,w)) +F = i [ ^l^ e ^/ 2 e^/ 2 tanh(^^)exp[zx(e^ -e^)] 
J An 2, 2 

x exp[-y(e pl + e p2 )] - -—. 

47TX 

Explicit evaluation gives 



o 



{e(w,w)} 



+F 



4nx y / x 2 + y 2 ' 



in agreement with ( 4.13|) . 



By duality, one gets similarly 



■<0|0> 



+ 



(4.12) 



«W,W)) + F = ]-[-- 7= 1 I ~ T— > ( I- 13) 

4tt yx x^x 1 + y 2 J 4ttx 

where we used the formula 

ds L ds^ i S_ 1 -S lelx{si _ S2)e _ y{si+S2) _ I . ,; I., |, 

a/sT y/s^ si + s 2 ^a; xa/x 2 + ?/ 2 

On the other hand, one easily gets from conformal field theory results that 

1 V 



^l°>+ =itmh ^. (4.15) 



The computation leading to ( |4.13| ) was done for y > 0. By considering y < instead, one 
also gets 

+ (0|0) F =-»tanh^—. (4.16) 
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In a similar way, we can consider the two point function of fermions. By considering 
the cases where yxyz > and yiy2 < 0, one also finds the result 

+ <AIA>, =4fl0tll AzA (4.17) 



+ <0|0} F 

The previous formula are compatible with the standard description |l| of mixed confor- 
mal invariant boundary conditions, which involves putting a "boundary conditions chang- 
ing operator" at y = 0. In the case from free to fixed spins, this operator has to be the <pi2 
operator of conformal weight yg. From this identification, it follows that ( [4.11| ) generalizes 
to 

+ ((3 2 n ■ ..0i\O)f _ -n TT f , A ~ Pj (Am , 

— ^M~ F — - 1 11 tanh —*—> (4 - 18) 

and that the scalar products coincide, formally, with the form factors of the (eg right 
moving) spin operator in the bulk. This ensures in particular that the partition fucntion 
of an Ising model with a region of fixed boundary conditions of length L inserted into a 
wall of free boundary conditions decays as L -1 / 8 . 



5. The case h a hb < 

As mentioned earlier, we expect the results of section 3 to hold for the case h a hb > 
only; physically, this corresponds to a situation where the boundary condition at x = — oo 
is independent of y. Instead, when h a hf, < and large enough magnetic fields, spins for 
y < and y > tend to be oriented in opposite directions, in particular they go to ±A4, 
where M. is the spontaneous magnetization, at x = — oo. There is thus a frustration line 
inserted in the system - in other words, following standard mappings [JTJ , a spin operator is 



inserted at x = y = 0. Of course, the spin operator on the boundary acquires a dimension 
1/2, so in fact we have a fermion line emitted at the transition. We thus expect that the 
problem of H — boundary conditions will be described by transition factors involving an 
odd number of particles; from that perspective, the ground state itself must be thought of 
as containing a particle of vanishing energy, and momentum equal to im (so the partition 
function of a model with a frustration line extending to x goes as e mx ) , ie a particle at 
rapidity 6 = which we denote by I in what follows. 

It is easy to check that this description is adequate in the massless case, and for 
boundary conditions a = +, b = — , ie the transition from "fixed up" to "fixed down". 
From conformal invariance, one expects @, 
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M>+ _H £ K*)>-+ = -^(i-4^). (5.1) 

Clearly, this cannot follow from formulas in section 3, since our scalar products F and G, 
which are functions only of squares of fields, vanish for h a = ±/i&. On the other hand, 
formula (|5.1| ) follows from the description with an odd number of particles. To see this, 
observe first that in the massless case, the particle at imaginary rapidity has vanishing 
energy and momentum, and merely keeps the correct parity of fermion numbers. Consider 
then the one point function of the energy. Two processes can contribute to it: in one case, 
the transition simply emits the particle at imaginary rapidity, that just "goes through" 
the energy insertion without being affected, so the process contributes (e)+ = — 4^7- The 
other process occurs when the transition emits a particle at real rapidity, which is then 
transformed into the particle at imaginary rapidity by the energy insertion. The required 
energy form factor is obtained by taking the usual limit: 

(Q\\e(w,w)\\f3,I) =u;e p/2 {exp [e?(-y + ix)] - exp [e (-y - ix)] } 

As for the scalar product +(/3|0)_, it should be of the form Xe^l 2 (A an unknown constant) 
since the frustration line emitted by the boundary corresponds to a fermion insertion, of 
dimension 1/2 (in other words, the partition function of an Ising model with a region of 
— boundary conditions of length L inside a domain of + boundary conditions decays as 
J d(3e@ exp [— e@L] = j;, as required from conformal invariance considerations JjJ). From 
this, we get the second contribution to the energy one point function 



/d(3 
— {exp [e^(— y + ix)] — exp [e^( — y — ix)] } 



i\u x 
tv x 2 + y 2 ' 

thus recovering ([5J]) provided A = — u. 

As a further check, we now compute the two point function of the energy with the 
H — boundary conditions. Conformal invariance gives the result 



(e(l)e(2))+_ = { ^L- (1 - 4 sin 2 0l ) (l - 4 sin 2 a 2 ) 

+ \ [1—4 sin 2 (ax - a 2 )] 5 — [l - 4 sin 2 (a>x + a 2 )] 1 , 

r z 7* ~\~ 4xix 2 J 



(5.2) 



where r 2 = (x\ — x 2 ) + (yi — y 2 ) , sin a = , x Several processes are now allowed, 

y ] /x 2 +y 2 

without or with real particles emitted from the transition. In the latter case, processes 
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where the real particles emitted at the transition are destroyed by only one of the energy 
operators are identical to processes appearing in the evaluation of (|5.1| ). The only new 
processes, which we call interacting processes, must add up to, using (|5.2|), 



1 



4tt 2 



4 sin 2 ai sin 2 «2 4 _._ 2 / \ , ' 



^ sin (an - ol 2 ) + —~ sin (ol\ + a 2 ) 



X1X2 r 2 r 2 + 4xiX2 

After a bit of algebra, this reads 

1 X1X2 (rf + rf) — 2x\X2y\y\ ~ 2x\x\ — 2xix 2 [x\y\ + x 2 yf) 
r 2 + 4xiX2 TT 2 r\r\r 2 



(5.3) 



where r\ = x\ + y\. On the other hand, there are two interaction processes. In the first 
one, the transition can emit a particle at rapidity fa, the first insertion of the energy 
can destroy the particle (3i and create one at rapidity fa, while the second insertion of 
the energy destroys the latter to replace it by the imaginary one. This process has the 
amplitude 

r dj^dfa^^ Pi _ _ e02 + 

J Z7T ZTT 

exp [e pl (-yi + ix±) - e p2 (-yi - ixi)] } 
{exp [e p2 (-y 2 - ix 2 )] - exp [e p2 (-y 2 + ix 2 )] } ■ 

After a bit of algebra, this reads 

If 2 yiy2-xix 2 -r\ 2 y x y 2 + x x x 2 - rf\ 



An 2 \r 2 + Ax\X2 rf r 2 r\ 

In the second process, the transition emits a particle at rapidity fa, the first energy insertion 
adds up a particle at rapidity fa and the particle at imaginary rapidity (the particle at 
rapidity fa just going through), and the two real particles are finally destroyed by the 
second energy insertion. This process has the amplitude 

r d^d^2 ePieP2 Pi _ zx2) + e02 ( _^ + 

J ZTT ZTT 

- exp [e pl (-j/2 + ix 2 )+ e p2 (-yi - ix 2 )] } 
{exp [-e p2 (-yi - ix{j\ - exp [-e^ 2 (-y x + ixi)] } ■ 

which is the same as the previous result, with 1 and 2 interchanged. Adding ( |5.4j ) and the 
similar result with 1 and 2 interchanged reproduces (|5.3| ) indeed. 
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We can now complete our study of scalar products by giving results equivalent to 
those of section 3 when h a h b < 0. Now only odd numbers of particles are emitted. The 
amplitude for an even number of real particles is 

b{ e 2n e,i\o) a = _ 

where the expression on the right hand side coincides with ( |3.30|) , while for an odd number 
of real particles we have the obvious generalization of (|3.30|) 



fa(02n+l, • ■ -dl\0) a _ T-r . K & + Z SUlh 6>i 1 



c n Y\_smh9if(9 



b (I\0) a 7111 " v w K6-isinhe iC osh(^ + f) 

* (5 6) 

ntanh — tanh — — — £ 
2 2 

The residue axioms fix c n again, except for an overall scale this time: c n = c (|) n . 

To fix this scale, we observe that at least the one particle form-factor should not 
vanish when h a = —h b , and that it should vanish when h a h b = 0, since then there is no 



frustration left. Recall that / itself goes as y7z| — /i 2 ,; we thus chose c n to be proportional 
to \fh a h b h / a 2 hb 2 ) i e we expect the one particle scalar product 

V h b~ h a 

b (9\0)a = ^v^ b ^^sinh9f(9 f b + lS ^ 9 Q , (5.7) 

2v/2 \fKi-hl ^ 6 -isinh^cosh(f + f )' 1 J 

where the field dependent prefactor and numerical factors are matched to the massless 
limit. Indeed, in that limit, h w e' 3 / 2 , giving rise to 

6</?|0>a = (/la " /lb) /3a) ^ _ ■ 

Finally, in the limit of large fields /i a , h b — > oo, h a /h b finite, this goes, using $(/3|/3&, /? a ) ~ 

2e i 7 r/4 e /3/2 e -^/4 e -/3 i) /4 j tQ 

& (/3|0) a = c^=^e^ 2 
Vh a hb 

In the case of opposite fields we recover the previous result for H — boundary conditions. 
Observe that, when h a = —h b , only the one particle scalar product is non zero, due to 



the terms y/n a — K b in the functions / ( since the only undetermined quantity is an overall 
normalization, we cannot change this without spoiling the results in the one particle case). 
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In the perturbative limit, when h a and /if, — ^ 0, one has 

sinh# 



b(9\0) a oc y/h a h b (h a - h b ) CQshdcQsh ^e _ W\ • ( 5 - 8 ) 

This can be compared with perturbative results. To do so, we have to slightly modify the 
action (|2.1|) to take into account the fermion line emitted at the origin. It is natural to 



represent it by a term ^h a hbi^>oa(0), where a is the same boundary fermionic degree of 



freedom as in ( |2.1|) , the term \/h a hb is dictated by a, b symmetry and dimensional analysis, 
and 1^0 is a fermion zero mode obeying (i/jqi^ l,r) ^ e mx . There are now two non trivial 
scalar products at second order: one does not involve insertion of this term, and thus is 
identical as the one we had for h a hb > 0, in agreement with ( |5.5| ). The second involves 
insertion of this term, plus a single integral along either the region y < or y > 0. This 
goes as (h a — h b )\/h a hb J_ dy((ijjL + ipR) (x = Q,y)ipL(x,y)), and the integral can be 
shown, using the propagators in the appendix, to agree with (|5.8|). 



6. Sinh-Gordon model. 

In this section we proceed to try to construct some of the transition factors for the 
sinh-Gordon model. Our purpose here is to gain control of the axiomatic relations satisfied 
by these factors in cases where the S matrix is non trivial, hence paving the way for further 
study of the sine-Gordon case, which is the most important for applications. As we will 
see, some complications arise when we wish to identify operators creating the need for a 
more systematic study. 

Let us first introduce some notation. The action is 



A 



dx 



dy 



1 



in 



t- cosh / 



+ A 



dy cosh - - 0o ) (x = 0,y). (6.1) 



The two particle S matrix and the reflection matrix for this model are well known. Setting 
sinh(| + ^) 



£(a) 



sinh(f -i^f) 

( r°° eft sinhf cosh(l - f)t 
Jo 



exp 



sinht 



sinh( T 



et 



in 



(6.2) 



< Re a < 4, 



they are given by 



R = 



£(fl)£(2-fl) 

^(1)^(2- 5/2)^(1 + 5/2) 
C(l - E(b))£{l + E{b))t{l - F(b))£(l + F(b)) 



(6.3) 
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1 

where B = 1+^/^ 1 an< ^ E,F are functions of b, also depending on the boundary 
parameters: they are given by E = and F = i^- in Ghoshal's notation [ pT2j] . 

In the following we will choose to work with the boundary action ( |6.1| ) having the 
phase 4>q = which leads to 77 = E = and simplifies the reflection matrix. 

Associated with the product £(1 + — 0) we define T(a) by 

cosh # , , & N 

/ 1 sin ( — ) 

t cosh 7T sinht Z7r 




with 9 



dt cosh ^ sinh | 
£ cosh 2 1 



Z7r — This function is constructed so that it satisfies 
T(a,6> + 27rz) = T(a,-0) ' 



(6.4) 



T(a,6>) 



(6.5) 



£(l + a)£(l-a) 

and has no poles nor zeroes in the physical strip. We will need this function in the next 
section to construct the form factors. Our choice of normalisation is such that 

T(a, 6>)Y(a, 9 + in) = sinh 9 + i cos (6-6) 



Boundary form-factors as well as transition factors can now be obtained by solving the 
set of axioms written in the previous sections. Let us start by investigating some general 
properties. 



6.1. Two particles sector. 

The two particle sector is where differences with the Ising model are emerging. Let us 
suppose that we have a set of boundary operators (local) denoted generically by k, then 
the transition form factor between states in theory a and theory b can be inferred to have 
the form 



,,-/.- ,n n , _ ,./.■ f if, . f ,n \ v ) F min{9l ~ 9 2 )F min (9 1 + 9 2 ) 
i 'ab\yi,V2) - C 2 Jab{9i)} ab [9 2 ) X 



Q\9 1 ,9 2 ). (6.7) 



cosh(^p) cos h( 



In this formula, F„ 



2 1 2 1 

is the usual function generating the sinh-Gordon S matrix 



Fmin(9) = A/^exp 



Mi 



B 



exp 



30 dt sinh(^) sinh [| (1 - f )] sinh f 
t sinh 2 1 

dt sinh(^) sinh [§ (l - f )] sinh \ 



sin 




sinh 2 1 



(6.1 
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All the dependence on the boundary is contained in the factors / a &. Q k is a polynomial 
we choose to be invariant under reflections, exchange and 2ni shifts. This leads to the 
following conditions on the f a b s 



fab(0) = R b (e)f ab (-e) 

f ab {9 + 2ttz) = R* a (0 + iir)Rt{9)f ah {9) 
for which a minimal solution can easily be found. Using the crossing 



(6.9) 



R i 



n 



9) = S(20)i?(i| +0) 



(6.10) 



we can write this last relation as 



In order to construct this function, we start with the functions u and v satisfying 

u(9)u(9 + in) 1 



(6.11) 



F rnin (29 + in) 
v(9 + 2m) = ^v(9). 

The function u will generate the S(29) and is easily found to be: 



(6.12) 



u(9) = exp2 



and v is given by 



fOO 



dt sinh ^ sinh [(l - f ) |] sinh f 



2 

sinh tcosht 



cos 



v(9) = 



T(F a ,9) 



(6.13) 



(6.14) 



?(F b ,ey 

The superscript refers to boundary conditions a, 6, corresponding to parameters F a , F b 
(since we have set 4>o = in the action, a and 6 correspond simply to different values of 
the coupling A). Now that the monodromy relations are satisfied, we need to multiply by 
other terms to get the correct reflection equation. This leads to the minimal solution 



Me) = 



sinh$ 



sinh 



9 , - 7T(1-F a ) 

2 "r 4 4 



sinh 



1 „- 7T(l+F a ) 
2 "r 4 4 



T(F a ,0) 
T(F b ,9) 



(6.15) 
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where we used the fact that u(6)/u(—6) = _ ^( 1 +- B /^( 2 B / 2 ) _ ^his minimal solution 
satisfies 



fab(e)fab(0 + i7T) = 



-4 sinh 2 9 



sinh 9 + i cos ZE |^- 



cosh(6>-z^)cosh(6> + ^) I sinh 6> + z cos ^ / F min (29 + in) 



-4 sinh 9 



(sinh 9 - i cos £§k) (sinh # + z cos 2f±) Fmm(20 + in) 

(6.16) 

This, added to the residue axiom leads to an equation for the polynomial Q k . Before 
embarking on the explicit form of these polynomials for specific operators, let us generalize 
these results to the N particle sector. 

6.2. N particle sector. 

Having understood the two particle sector the general structure emerges nicely and 
turns out to be a very simple generalisation of the two particle case: we have 



i=i 



><n 



F m in(,@i 0j)F m i n i y 9i Oj J I 



(6.17) 



i<j 



cosh(^)cosh(^ i ) 



Q (#1, • • • , n ). 



By construction, this satisfies all the exchange and reflection relations. Here only the 
residue equation needs to be fixed: it leads to an equation for the polynomial Q k in the 
form 

Q k (0 + Z7T, 0, 3 , • • • , n ) = P(0\0 3 , ■ ■ ■ n )Q k (0 3 , • • • , 0„), 



(6.18) 



with P given by 

P(0\0 3 ---0 n ) 



1 



sinh 



7rF a \ ( ixF h 
sinh 9 — i cos | | sinh 9 + i cos — — I x 



n 

i=3 



sinh(# — Oj) + sinh 



2 

inB 



sinh6> + i cos J ^sinh6> 



n 

i=3 



sinh(6> — Oj) — sinh 



inB 



sinh(6> + Oj) + sinh 

. 7TF b \ 
l cos X 

2 J 

sinh(6> + Oj) — sinh 



inB 



(6.19) 



inB 



and we have adjusted the constant c n in order to absorb the remaining constant factors. 
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6.3. Operator identification. 



These equations are somewhat similar to those found for operators in the bulk |T3 
There are some essential differences though, which make the operator identification more 
difficult. For instance, in the bulk, an argument based on Lorentz invariance can be made 
to determine the degree of the polynomials Q fc , depending on the spin of the operators: 
it is not clear how to extend this in the boundary case. In order to do the operator 
identification, we will use below the fact that one can analytically continue the form factors 
in the coupling constant B —> l +—/3 , leading to the so-called roaming trajectories between 



minimal models JI^]. Then, as (3q — ► oo, at least for some simple operators, we should 
recover the form factors of the Ising model obtained in the previous sections. As we will 
see, this is not sufficient, in fact it allows to determine only part of the polynomials Q k . 
Another way can be devised to identify operators by taking the massless limit: if we know 
the boundary dimension of the operator we can identify formally the form factor series 
with that of the bulk. This will allow to determine the maximal degree of the polynomials. 

A specific property of the boundary case is the following: requiring the polynomials 
to be symmetric under exchange, and invariant under reflection, leads to the fact that they 
should be symmetric polynomials in the variables x% = cosh^. Naively, one could then 
be tempted to construct homogeneous symmetric polynomials in these variables, but this 
would be wrong: since the polynomial has to be symmetric in e 9 and e~ 9 we can form all 
invariant polynomials in these two variables. The sum is cosh# but the product, which 
is another symmetric polynomial in these variables, is 1. This means that we should not 
expect the polynomials to be homogeneous. 

We will not carry this study much further here, and contend ourselves by giving a few 
examples, mostly in the case of scalar products, ie transition form-factors for the identity 
operator. Since we choose our form factors to be normalised to 1 in the zero particle sector 
(we divide by the expectation value) , we obviously have that the zero particle form factor 
in simply 1. For the two particle scalar product, we have the equation 

Q Id (9 + m,9) = 1. (6.20) 



where we adjusted the constant c 1 ^ such that 

irF, 



Id COS 2 - COS 



(^) 



(6.21) 



24 



Obviously there are more than one way to satisfy this equation. Choosing Q Id (9i, 62) = 1 
would seem, at first sight, to be the correct solution for the scalar products. Following 
let us take the roaming limit B = 1 + -P with f3 Q very large. Then we observe that 

F min (9)^-ie-^/ 2 smh 6 - 

2 (6.22) 

^yinh-G falsing _ 

Thus the scalar product becomes that of the Ising model provided we choose iF a = —9 a . 
On the other hand, if we take the massless limit and choose boundary conditions a and 
b to be of Neumann and Dirichlet type respectively, then from conformal field theory the 
situation is described by the insertion of an operator of dimension h = 1/8. Then we can 
compare our choice for Q Id with the expansion of massless operators in the bulk having 
h = 1/8. The conclusion, at this 2 particle level, is that there should also be a contribution 
proportional to o\ = J2i=i x i = Xli=i cos h^- This term would not modify the residue 
equation. Thus, based on this analysis we expect 

Q Id (xi,x 2 ) = Ao\ + 1. (6.23) 

The main difficulty is that we have no means to fix the relative normalisation A and we 
are left with a free parameter. 

The four particle contribution is constructed using the residue equation, we have 

Q Id {-x, x, x 3 , x 4 ) = P(x\x 3 , x 4 )Q 7d (x 3 , x 4 ) (6.24) 

which allows to construct the four particle scalar product using the property of symmetric 
polynomials § 

a£ +2 (a;, -x, xi, x n ) = a% - x 2 a%_ 2 . (6.25) 
The part coming from the "1" can be fixed uniquely again 

Q Id,1 {x U X 4 ) = 4 7 ((73(7 2 - 7 2 C7 3 ) + 4~f(K a K b - l)(-£7 2 0"i + 7V1) 

(6.26) 

- 2i(K a - Kh ){al - a 3 a 1 - 4(1 + 7 2 )a 4 - - 2 7 2 a 2 + 7 4 ). 

Here X{ = cosh^j, 7 = sinh ^5 anc [ _ cos ZL^k/i. "We observe that as expected, this 
part of the solution also collapses to that of the Ising model in the roaming limit. 



6 Here we use the notation for the elementary symmetric polynomials defined by the gener- 
ating function fj fe (l +txk) = 'Yl lk t k <7k- The superscript indicates the number of indeterminates. 
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For the other part, there is an ambiguity for the term of degree 6 coming from the 
kernel of the residue equation 



Residue (6, 9 + ztt, 3 , 4 ) = (6.27) 

and in general we will need physical arguments to choose the correct solution. We believe 
a full study of the boundary operator content of this integrable theory could lead to a 
complete solution but felt it is beyond the scope of this paper. 

Finally, we dicuss the boundary form-factor for homogeneous boundary conditions 
(a = 6), in the simple case of the trace of the stress energy tensor. The simplest guess for 
two particles is 

Q e (d x , d 2 ) = cosh 9 ± + cosh 9 2 

= 2 cosh( ) cosh( j. 

Here the addition of a constant term would spoil the residue property needed. The normal- 
isation is fixed by taking the roaming limit, where the form factor should be identical with 
the one of the Ising model (computed explicitely using the mode decomposition given in 
appendix A). Moreover, the massless limit reproduces the bulk form factor up to reflection 
matrices. 

The four particle form factor is determined using the recurrence relation of the sym- 
metric polynomials (|6.18|) and we find 



Q e (6 1 , ■ ■ ■ ,0 4 ) = [4 7 («a " 1)(-*2<ti + 7V1) + 4 7 (a 2 a 3 - 7^3)]^. (6.29) 

It is clear that in the roaming limit, this form factor vanishes, as expected since in the 
Ising model the four particle form factor of the energy-momentum tensor is zero. Moreover, 
assuming that o\ must be factorised, the ambiguity in this polynomial can be fixed since 
by taking the massless limit we know it must have at most degree 6. 

We believe that these few examples are sufficient to show that the method is applicable 
to any integrable model, and in fact is very similar to the usual problem of finding form 
factors in the bulk. However, there is a need for a systematic way to identify operators 
at the boundary and determine the free parameters discussed above. We now conclude by 
some applications. 

26 



7. P(t) in dissipative quantum mechanics 

The applications of our formalism to 2D statistical mechanics problems are somewhat 
obvious. We rather concentrate on quantum impurity problems, and discuss now a stan- 
dard problem of dissipative quantum mechanics, the computation of "P(t)" in the double 
well problem. In the so called ohmic regime, the two state problem with dissipation can 
be mapped on a single channel Kondo model ||, with hamiltonian 

ro i , , 

H x = dx-[Il 2 + (d x (j)) 2 ] + \5(x)(S + e l ' / ^' i, + S-e- lV ^^). (7.1) 

J-oo 2 

Here, S± are spin 1/2 operators, the values up and down corresponding to the two states 
of the system. The dissipation is characterized by a dimensionless number, which can be 
taken to coincide with the conformal weight h = g of the boundary operator. 

The physical quantity of interest, P(t) is defined as follows: assuming that the spin 
has been fixed in the up state up to time t = 0, P(t) is its average value as a function of 
time after turning on the dissipation. Mathematically, introducing the ground state of the 
theory where the spin is not coupled to the heat bath, |0)a=o, we need to evaluate the one 
point function 

P(t) = {{l\S s \Sl), (7.2) 

where |fi) = e~ zHt \0), |0) is the tensor product |0) a=o ® |+) and H is the hamiltonian ( |?~T] ) 
with dissipation. Equivalently, we can write 

P(t) = (0\e lHt S z e- lHt \0). (7.3) 

We denote the quantity inside the bracket as S z (t) (evolution with respect to the hamil- 
tonian ( |7.1|) is implied). 

We adress the problem in the formalism described previously. Introducing a complete 
set of eigenstates of H\ on the left and right hand side of flO]), the evaluation of P{t) 
requires the knowledge of the scalar products of these eigenstates with |0)a=o <8> |+ >• To 
diagonalize the hamiltonian, we then use integrability of the anisotropic Kondo model, and 
the resulting massless scattering description, involving solitons and antisolitons (of charge 
±1) and breathers. 

There is a slight ambiguity in this approach, because the scattering description is 
based on an infrared picture, where the spin is always screened, for A > 0. For any A > 0, 
the spin can be "extracted" from the asymptotic states: for instance, the spin correlators 
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follow simply from the correlators of d x (f>, which can be computed using the form- factors, 
without any reference to the spin. At A = 0, it is not clear "where" the spin exactly is in 
our description. Nevertheless, one can write axioms for the scalar products of multiparticle 
states at A ^ with |0)a=o <8> |+ >, and of course they are of the same type as the ones 
written before for any pair of couplings A, A', so it seems our approach should work in that 
case too. The effect of the spin |+) in |0) is now simply taken into account by considering 
that |0) has a charge equal to one, ie the only non vanishing scalar products will be of the 
type ;r +1 - ei (/Wi ■ - -/?i|0>, with £ = 1. 

We discuss first the case g = 1/2, where the Kondo problem essentially decomposes 
into two decoupled Ising problems, with boundary field h oc A JT^] @. At this value of g, 
the spectrum is made of a soliton and an antisoliton, and the R matrix is 

R + _ = R-+ = i^^, R++ = R- = 0. (7.4) 

and only the two particle form-factor of the current is non zero 

(0\d x <f>(x = 0,t)||/3i/3 2 > ei£2 oc 5 £l+e2 e ie ^/V 2 / 2 
[1 + R + -((3 1 )R+-((3 2 )} exp [-it{e& + e p2 



(7.5) 



To proceed, we have to compute the one point function of d x <p by inserting complete 
sets of states. Many processes will contribute. In the simplest, a soliton is emitted at the 
transition, acted on by d x (p that transforms it into another soliton, and this second soliton 
is destroyed at the second transition. The amplitude for this process is 



1W 



<0|/3 2 )+ +(^|0> +</W(Wi>J- (7-6) 



This process gets "convoluted" with many others, where an additional arbitrary even 
number of particles is emitted at one transition, simply "goes through" , and is destroyed 
at the other transition. This means one has to replace (0|/?2)x \ i n ( |7.6| ) by 



(o\fh)t nm+ 1 ^f i <oi/w 2 >r + r + </wiio>+.... (7.7) 



7 In the limit A — > oo, the change of boundary conditions corresponds to the insertion of 
an operator with dimension h = |, twice the dimension of the Ising spin (for arbitrary g, this 
dimension is f ; in particular it is 4 at the isotropic Kondo point [17] 



28 



The manipulation of this sort of terms follows from the general principles used in section 
3.3: we simplify them by moving contours and using closure relations. Here, we move 
the Pi, fa integrals to Im Pi = —in and Im (3 2 = +in; by doing so, we do encounter 
singularities of the form-factor of d x (p, and kinematical poles. Forgetting these singularities 
for a while, crossing produces a 6 (Pi — P2), which gives a vanishing contribution because 
(P\d x (p\P) = 0. The kinematical poles, together with the coincident terms in crossing, 
give rise to terms like S(j3i — (3%). By proceeding similarly with the new integrals, we can 
succesively eliminate all the integrals, and we are left with the contribution of the poles of 
the form- factor of d x (j). Going back to ( |7.6| ) and ( [T.7|) , when we move the contours of Pi 
and fa integrations, there is a pole for Pi = Pb — ^ and fa = Pb + -y- For each of these 
poles, we are left with some integrals which we eliminate in the same fashion as before: all 
what remains at the end is proportional to the residue of the d x <p form factor at the poles, 
ie a term e~ 2Tbt , up to some numerical, time independent constant. 

There are two other types of terms contributing to the one point fucntion of d x <p. 
In the first type, the first transition emits two solitons and an antisoliton, a pair soliton 
antisoliton is destroyed by the d x (p insertion, and the other soliton by the second transition. 
In the second type, the d x <p insertion emits a pair soliton antisoliton. These processes 
get convoluted with many others as before. By moving the integrals, we find again a 
contribution proportional to e~ 2Tbt . 



Using now the relation between d x cj) and the spin |T8 



S z (t)~S z (0)= I d x cf>(x = 0,t')dt', 
Jo 



we find immediately 

P(t) = e- 2T »\ g=\, (7.8) 



a well known result |T9[ . 

Remarkably, the only contribution to P(t) comes from the poles of the R matrices, 
or equivalently of the non kinematical poles of the transition factors. Since the contour 
manipulations did not involve the explicit form of G (nor the proof that two different ways 
of computing the form- factors give the same result, see section (3.3)), it is very likely that 
this feature generalizes to the case g ^ h. 

Assuming thus that P(t) is completely determined by the poles of the transition fac- 
tors, that is the poles of the R matrices, we can now write a very natural set of conjectures 
for its general form. 
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First, recall that for g > h, the sine-Gordon spectrum has no bound states. The R 



matrix for solitons antisolitons is independent of g: R± 
We thus expect: 



i tanh 



P-Pb 



m = E 



n=l 



a n e- 2nT »\ gy 1 - 



0. 



(7.9) 



with ^ a n = 1 . There are an infinite number of terms here compared to the g = \ case 
because now all the (U(l) neutral) In particle form-factors are non zero. 

On the other hand, for g < », there are bound states in the spectrum. This means 
that the transition factors and the R matrices now have non trivial poles. According 
to our conjecture, these poles contribute to P(i) by values of having a non zero real 
and imaginary part, and we thus expect that P(t) will develop oscillatory components on 
this side of g, in agreement with known numerical results, and an analytical argument in 
perturbation around g ||20|| . At large times, we expect the dominant contribution to come, 



as usual, from the one-breather term. Recall the one breather reflection matrix [HB 

tanh 



R(J3) 



iirmg 
4(1-3) 



tanh 



+ 



(7.10) 



together with fii 
mass. 



2 sin 



2(1-9) 



4(1-3), 

the ratio of the breather mass parameter to the soliton 



This should give the leading behaviour of P(t) (there is no factor 2 in the exponent 
now because a single breather can be destroyed by d x <p) 



P(t) oc exp 
In particular, setting g 



-2m sin 2 



Tig 



cos 



2(1-0) 

1 — e one gets 

P(t) oc e- 2Tbt cos {4T b net) 



tTh sin 



Tig 



(1-9)1 



(7.11) 



(7.12) 



The relation between T b and the bare parameter A is Tt oc A 1 / 1-9 - the coefficient of 
proportionality can be found inpTJ. Expression ( |7.12| ) agrees with the expansion near 
g = 2 in ||20|| . The expression (|7.11| ) is new as far as we know: in standard notations, we 



predict the ratio of the period of oscillations to the damping factor to be 

O ng 

— = COt —. r . 

r 2(1 -g) 



(7.13) 



It would be interesting to test this numerically (in fact, ( |7.13j ) agrees very well with mea- 
surments that appeared recently in a paper of K. Voelker |p2| ). 
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8. Conclusion 



It is clear that more work is necessary to fully determine the transition form-factors, in 
particular the boundary form-factors as well as the scalar products, for general integrable 
theories. This seems to be a purely technical matter, and we believe that all new qualitative 
aspects are well described by the examples discussed here. 

To conclude, we would like to make some general remarks, illustrated in the Ising case. 
From the knowledge of the transition factors, we have in fact obtained an expression for 
any state of the theory with boundary condition a in terms of the states with boundary 
condition b. In particular, we have 



|0>c 



6<oi°>«e /n^ G( ^ w ^ 

n=0 i=l ^ ' 



Using that ground states are normalized, the scalar product of the ground states follows: 



1 =E/n§(4i«(*-.. -.*)! 

n=0 J i=l v ' 



»<0|0) a | 



2 



It is not clear how useful this general expression is. For instance, in the massless case and 
for free and fixed boundary conditions, this is nothing but the (R part) of the spin spin 
correlator at coincident points in the bulk, a manifestly divergent quantity. 

It is also interesting to get back to the direct channel, and express the boundary 
state in the case with two different boundary conditions. By matching the expressions of 
correlators in both channels, one finds the expression 



B) ba = exp || ~ [K a {6) + K b (6)] Z*(-6)Z*(9) 



|0>, 



where the poles at opposite rapidities 6i = —9j are regulated by taking the principal part. 
Going to the crossed channel, the contours have first to be completed by turning around 
the singularities from above or below, the residues completing the term in the exponential 
to give either K a or K b . The contours can then be moved by ±^ to give integrals in terms 
of F or G, reproducing the two possible expressions for correlators as discussed in details 
for the energy in section 3. 
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Appendix A. Perturbative computations. 

In this appendix we show some of the direct checks done on the form factors discussed 
earlier. If we use the action of the introduction, we can decompose the fermions as 



4> 



d9[cue e/2 a(9)e imx sinh 0+imt cosh e + we^V {9)e 



imx sinh 6 — imt cosh 9l 



-oo 
oo 



d6[ue- e/2 a(8)e irnxsinh9+imtcoshd + uje~ el2 a\9)e 



(A.l) 



imx sinh — imt cosh 9] 



where we have used uj = e l7r / 4 . We evaluate the scalar products between different states 
in the simple case where the fields h a = and hb are small. The states at t = — oo are 
in the Hilbert space of the theory a, which is at zero field, and those at t = +00 are in 
the Hilbert space of theory b. Then from standard perturbative arguments we obtain to 
second order (which is the first non trivial one) 

7 2 roc i>oo 

b (eie 2 \\o) a <x-± / d yi dy 2S gii( yi - y 2 ) / de[de' 2 e- m yi cosh ^~ my * cosh 

z Jo Jo 



2 cosh 



n 

2 ~ l A 



+ 2i2S(^)oosh[-^+< 4 



n 



. 7T 



2 C osh f ^ - i-j + 2R* (9[) cosh [f+i^ 
x (0\a(9 1 )a(9 2 )a\9[)a^9' 2 )\0) 



0'i 



(A.2) 



7T 



where Rq is the reflection matrix at zero field. The integral (|A.2| ) are easily computed by 
using the mode decomposition for the fermions: 



b (M 2 iio) a0 c^ n sini1 ^ 



'~~ 2 cosh^ cosh(-^- — i- 



tanh 



9\ — 9 2 



tanh 



0i +9 2 



(A.3) 



We can do a similar computation for the four particle scalar product to fourth order and 
we find 



,(9, ■■ -8,110)0 oc hif[ Sinl1 "' 



4 cosh6*i cosh(-^- — i 



Yl tanh 



A' i<j 



9i — 9j 



tanh 



9i + 9j 



• (A.4) 



34 



